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Abstract The Indirect Detection of neutralino Dark Matter is most promising through annihilation channels 
producing a hard energy spectrum for the detected particles, such as the neutralino annihilation into 
Zh. A cancellation however makes this particular annihilation channel generically subdominant in the 
huge parameter space of supersymmetric models. This cancellation requires non-trivial relations between 
neutralino mixings and masses, which we derive from gauge independence and unitarity of the MSSM. To 
show how the cancellation overshoots leaving only a subdominant result, we use a perturbative expansion 
in powers of the electroweak/supersymmetry breaking ratio mz /m^. 
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1 Introduction and motivations 



There is no doubt that standard models of particle physics 
and cosmology alone cannot describe the full wealth of 
observational data recently collected on a wide variety of 
length scales. Ad-hoc as it may seem, the Dark Matter 
(DM) hypothesis [1,2] is probably part of the minimal set 
of extra ingredients needed to account for the increased 
gravitational self-attraction of matter on scales ranging 
from galaxies to the full visible universe. More exotic in- 
gredients like repulsive Dark Energy, or modifications of 
gravity itself, might also become necessary to cope with 
the apparent acceleration of the universe. In the absence 
of a convincing unified theoretical solution to these both 
issues, experimental searches are the only way to prove 
the validity of hypotheses like the existence of a DM par- 
ticle. For instance, the first issue would be settled if a new 
particle were found and its non-gravitational interactions 
measured to be compatible with the Cold Dark Matter 
relic density required by Cosmic Microwave Background 
and large scale structure formation [3,4]. However, these 
interactions are then by definition weak, and we should 
not be surprised that their evidences are extremely hard 
to obtain, much like for Pauli's neutrino. Debates like the 
one around DAMA's claim [5-12] for direct detection of 
DM are illustrative of this difficulty. This is why it is cru- 
cial to be able to cross-check and understand results in as 
many different ways as possible, for which a definite and 
well motivated DM theoretical framework is necessary. In 
this work, we shall keep with the well-studied supersym- 
metric lightest neutralino (x) of mSugra or MSSM models. 



A particularly crucial cross-check would be the Indirect 
Detection of the neutralino annihilation products, which 
would show that such annihilation did indeed occur in the 
past, froze out at some point and re-started in hot spots 
like the galactic center or the solar core, where dark mat- 
ter later accumulated. However, to identify an Indirect 
DM signal and possibly determine the neutralino mass by 
looking at fiuxes of e.g. photons or neutrinos from such 
hot spots, it is essential to be able to distinguish that 
signal from the standard but poorly known astrophysical 
background. These being characterized by energy spectra 
with fairly universal power laws, Indirect Detection will 
be most successful when neutralino annihilation proceeds 
through primary channels which provide secondary pho- 
tons or neutrinos with the hardest possible spectra, and 
a sharp energy cut-off around the neutralino mass. From 
fig-ffl (discussed in Appendix [S]) , the most promising an- 
nihilation channels are into a r+r~ pair, into two gauge 
bosons (xx ~* W^W~ , or ZZ which has the same shape) 
or into one gauge boson and a Higgs-Englert-Brout (HEB) 
boson [13,14] ixx^Zh). 

However, these fairly universal spectra need to be weigh- 
ted by the actual model-dependent branching ratios to 
give the final indirect DM detection signal. It was noted 
long ago [15] that the Zh channel is then suppressed, 
which can be numerically checkecQ using the DarkSusy 
(3.14.02 version) [16] and the Suspect (2.003 version) code 
[17]: the top plot of fig. [2] typically shows a suppression 
by three orders of magnitude for various mSugra mod- 
els with tan/3 = 10. To qualify this suppression, let us 
start by contributions we expect to be dominant, namely 
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Figure 1. The shape of differential neutrino (top) and gamma 
(bottom) fluxes coming from neutralino annihilations into 
various primary channels. The normalisations (depending on 
branching ratios) have been arbitrarily rescaled to compare 
the shape of all channels; the dependence on (1 TeV here) 
is weak. 



those from SM particle exchanges (in our case an s-channel 
Z) since superpartners are necessarily heavier. As seen in 
the bottom fig. [2l this contribution not only dominates, 
but even overwhelms the total cross-section. We there- 
fore need to understand a double suppression, that first 
cancels this large contribution and second brings the total 
Zh annihilation below other channels. The cancelling con- 
tribution necessarily involves non-SM particle exchanges, 
which seems contradictory with the fact that on fig. El the 
cancellation gets better with increasing mg and mi/2, i-e- 
for maximally broken supersymmetry. 

To be as general as possible, this cancellation should 
be checked in a supersymmetry breaking independent way. 
This is done in fig. [3] in the more general MSSM, while 
keeping the usual GUT relation Mi = | tan 6*1^-^/2 — 
O.5M2. A cancellation up to three orders of magnitude 
thus seems a generic property of every broken supersym- 
metry theory. 

To be more concrete, let us now focus on the particular 
mSugra model with mp = 3000 GeV , mi/2 = 800 GeV , 

= 0, tan (/3) = 10 and ^ > 0, marked by the black star 
in fig. [21 The annihilation cross-sections at rest are: 
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Figure 2. Neutralino annihilation branching ratios to Zh for 
various mo and mi/2 values: Z-exchange only (bottom), Z and 
X exchanges (top) 



/i (GeV) 




150 200 300 400 500 600 700 800 900 Mi 

{GeV) 

Figure 3. Contour plot of in the MSSM; white 

dots correspond to mSugra models of fig. [2] 
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The annihilation is dominated by the tt channel, 3 or- 
ders of magnitude larger than the Zh one [18]. However, 
when restricting to the Z exchange diagram, they are com- 
parable, with Zh slightly larger. This is easily understood 
by exhibiting the couplings and kinematic factors in the 
amplitudes: 



h 




[XX 



tt] 



cos-^ 6w 



X 



A ixx ^ Zh) 2 cx 



where the XX^ coupling Ofj^ is defined in terms of the 
neutralino mixing matrijo N (see Appendix ID.ip : 



(1) 



It remains to understand why other exchanges can- 
cel the Zh channel and not the ti one. For the latter, 
the t-channel sfermion exchange can be made arbitrarily 
small by taking large enough mo, so that the main contri- 
bution proceeds via a (SM) Z boson exchange, conform- 
ing to naive expectations. However, two other diagram 
are involved in the Zh annihilation channel: (1) an s- 
channel pseudoscalar A exchange, which can be neglected 
for large mo (in what follows we will always work in this 
pseudoscalar decoupling Hmit l(47| ). and (2) i-channel ex- 
changes of the 4 neutralinos, which cannot be decoupled 
because of strong links between the couplings involved in 
each diagram. These links appear in the expression of the 
annihilation amplitudes [19, 20] derived in the next sec- 
tion: 



A{xX^Zh),^-^^P,,xOf, 



A ixx ^ Zh)^ 



X 

cos^ t/w rnj^z 
ig^^ ml 
cos^ Ow fn^z 



(2) 
(3) 



^ 2m?, + 2m?,. — m? — m?r, 

where the XiXjh coupHng O^j is defined in Appendix l(i9| 
as (where the decoupling condition fTZl) has been used): 



icwN,2 - swNa) {spN-ji - c^iVjs) + (« ^ j) (4) 



It is clear that the second amplitude ^ cannot easily be 
neglected and might turn out comparable with the first 
one 1(2]). However, it is less clear why both should cancel 



^ In what follow, we will always work with the hypothesis of 
CP violating phases absence, in order the N matrix to be real 



with high precision, especially given as different-looking 
coupHngs as (U) and 

A toy example of the above cancellation is provided 
by the annihilation of a spin singlet ti pair into Zh in a 
Standard Model without SU{2). One may first be puzzled 
to get a cancellation between an s-channel Z-exchange, 
which only involves gauge couplings, and a t-channel top 
exchange, which involves an a priori independent Yukawa 
coupling. However, one soon realizes that the existence of 
the s-channel requires both spontaneous breaking of the 
gauge symmetry for the ZZh vertex, and an axial coupling 
for the tiZ vertex. This last coupHng excludes contribu- 
tions to the top mass other than the gauge breaking y{h) 
one. The two channels q-^qmz oc q—y are then both 

proportional to gjiji) and may cancel. In contrast with 
this simple case, the neutralino annihilation studied be- 
low is complicated by the presence of another source of 
mass, namely the SUSY breaking Majorana mass terms 
for the gauginos. 

To analyse this cancellation, we start in section 2 by 
deriving the relevant amplitudes. In section 3, gauge inde- 
pendence is used to draw a first link between the couplings, 
which is shown to follow from the gauge invariance of the 
mass matrix. Unitarity at high energy is then used in sec- 
tion 4 to derive a second relation, which is combined with 
the first one to show that a cancellation is possible. In 
section 5, the structure of Rayleigh-Schroedinger pertur- 
bation theory is then used to show that this cancellation 
is stronger than expected. 



2 Amplitudes for Neutralino Annihilation at 
Rest 

Let us start by deriving the polarized amplitude A (Zh)^ 
for neutralino annihilation into Zh via Z exchange. By 
construction, neutralinos are their own antiparticles, so 
that an initial pair of lightest neutralinos at rest is nec- 
essarily in an antisymmetric spin singlet state. The final 
state containing a HEB scalar, the outgoing Z boson po- 
larization needs to be longitudinal, and can be chosen as 
z-axis. It is then more convenient to use helicity ampli- 
tudes [15, 19] than unpolarized cross-sections [20]. The 
Feynman diagram and rules defined in the Appendix [D] 
give the amplitude 



X,P 



h,k' 



XiP 




A{xx^Zh)z 



(see Appendix El 



2^"? (x(p)7^5Xb)) (5) 

q 
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where x = w is the incoming neutrahno, x — v the outgo- 
ing one in an arbitrary choice of arrows directions, u and 
V are external Dirac spinors at rest in the chiral basis (see 
Appendix [BJ , namely: 



or = 



in agreement with ^ and existing results [19,20]. 

This amplitude (flOl) has to be compared with the result 
of a similar computation for the annihilation into ti: 



cos^ tiw rni 



(11) 



for a spin up or down along the z-axis, and similarly for 
the "antiparticle": 



Thanks to the Majorana condition (|40| . descriptions of 
the same neutralino as a particle with or as an "an- 
tiparticle" with rjsi are equivalent provided: 



-ia £,s 



(6) 



The polarization of the longitudinal Z boson being e'^{k) ~ 
(fc^, 0, 0, k^)/mz, and the initial momentum at rest simply 
= (2m^,0), the polarized amplitude ([5]) is: 



where f3t, 



1 — m^/rri^ and r3=l is the weak isospin of 



the top quark. Notice the diflerent power of m^, favouring 
the Zh channel for large masses. 

We have seen in the introduction that t-channel neu- 
tralinos exchanges should reverse this conclusion. Follow- 
ing the same path, their contribution is 



X,P 




X,P 



Z,k 



■Ki+'^fi) (7) 



(fc^O,0,A:°) 



4 



(8) 



Aixx^Zh) 



li 



X (Xb')(^ + m^j7"75xb))ep {k) 



with Q'' = p'' — /c^. After some algebra on the <-channel, 
and adding the u-contribution (same diagram as before 
with p and p' interchanged), one finds: 



Notice the time-like structure of the initial state vector 
(u^i = (2, 0, 0, 0): a purely axial coupling talks only 

with the scalar part of the two spins at rest. Notice also 
the disappearing of the Z pole at = mz /2. 

Expressing k^ = m-^^zh in terms of the conventional 
kinematic factor 



P. 



Zh 



{mh + mzY 
4rn^ 



mz) 



4m2 



1 



» mz 



and using the definitions l(43|) and l|46|) of couplings in 
terms of neutrahno mixings given in the Appendix [D] we 
finally find: 



cos^ 



w m% 



-P. 



Zh 



(9) 



To get the amplitude with reversed hehcities, we just 
need to replace (1,1) components of the PauH matrices 
in |[7]) by (2,2) components, and take the extra sign from 
the Majorana condition ^ into account. Thanks to this 
sign, only an antisymmetric initial state can contribute, 
and gives with the correct normalization factor: 



Aixx-^Zh), 



cos^ 



2 PzhOii 



(10) 



Aixx^Zh)^ 



'■V2 



^9 V xo ^ 

cos"^ Ow m^z 



(12) 



2™^ 



The problem is now to understand how the amplitudes 
(fT2l) and (fTO|) cancel with the 10^'^ precision shown in 
the introduction. This can happen only if the sum of sec- 
ond lines of (fT2|) which contain four powers of the neu- 
tralino mixing matrix N somehow reduces to Of^ ~ N'^ 
as a consequence of some symmetry. We already noticed 
that supersymmetry had to be maximally broken for the 
cancellation to take place. We are thus left with gauge 
invariance which is investigated in the next section. 



3 Gauge Independence and Gauge 
Invariance of the Mass Matrix 

The previous computations were performed in the unitary 
gauge, where massive gauge fields have completely "eaten" 
a Goldstone boson. One way to obtain non-trivial relations 
of the kind we seek is to work in the i?^-gauge family of 't 
Hooft [21], and require independence of the result on the 
gauge-fixing parameter ^. 

Let us first notice that the neutralino exchange dia- 
grams are gauge independent by themselves. Indeed, in 
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a vector supermultiplet, only the bosonic gauge field is 
gauge dependent, and not the associated gaugino. More- 
over, higgsinos are associated with the real part of complex 
scalars, which are also gauge independent. 

We can thus concentrate on the Z exchange diagram. 
When going from unitary to i?5-gauges, this contribution 
splits in two ^-dependent diagrams: one with Goldstone 
boson exchange, and another with the Z-boson exchange. 
To exhibit the cancellation of their ^ dependence, it is 
convenient in the Z propagator 



_ (1 - Ci qtiQp 

to decompose the longitudinal part as: 
m|(l-0 _ 1 



1 



m|) (g2 



ml) 



(q2 - 



(13) 



The first term is nothing but the longitudinal unitary 
gauge propagator, and the second exhibits a fake pole at 
the Goldstone mass = S,m\, with a wrong sign [22]. 
The ^-dependent Z exchange Az(^ = Az + Azg corre- 
spondingly decomposes into the previously obtained Az 
^ and a Goldstone-like amplitude with gauge couplings: 



A (xx ^ Zh) 



.C^C 



z /_ 



ZG 



ix {p)4i5x{p)) 



q-e{k) 
'2 - ^to| 



which, using the mass-shell condition for the initial state 
IbX = 2toxX75X, becomes: 



A {xx^ ^MzG = 



2 cos^ Qy^ mz 



We now turn to the genuine Goldstone boson exchange. 
The ZZh coupHng is replaced by the GZh one: 



G 



C^i^ = iCa {k - k') where Cg = 9/2 cos9w (15) 




and the scalar propagator is simply «/(g^ — C"^|)- But 
things are more subtle for the Gxx vertex: even if the neu- 
tral Goldstone boson is part of the Z boson in the unitary 
gauge, it does not have the same coupling to neutralinos: 



G 



XiXi 



2 cos 6* 



w 



(16) 



with 



For the Hghtest neutralino annihilation, i ^ j ^ \, and 
since the coupHng is imaginary, the Goldstone only couples 
to the axial part of the neutralino. The amplitude for the 
Goldstone exchange diagram is thus: 



A ixx ^ Zh)^ = -C^^Cg X {p') 15X ip) (18) 



{q + k') 



g2 - 



e^'{k) 



Recalling kinematics {q = k + k'), the polarization con- 
dition {k.e{k) = 0), and couplings definitions, we finally 
have: 



Aixx^Zh), 



ih;^ X15X 



q-e{k) 
g2 - 



(19) 



Now, comparing (|T4l) and l|19p . gauge independence re- 
quires a relation between couplings: 



^11 , 



,^11 



mz 



which by |(T|) and (fT7|) can be expressed in terms of the 
neutralino mixings and masses: 



14 



mz 



(20) 



- {N12CW - Niisw) (spNii + C13N13) 



This relation can be extended for the annihilation of an 
arbitrary pair of neutralinos Xi ~ Xjin the same channel: 



Of 



mz 



which is a short version of the rather non-trivial identity: 



"Xi 



mz 



(21) 



Although not completely identical, this relation bears 
similarities with the combination of masses and mixings 
involved in |(3]). It is therefore interesting to notice that 
it only involves the neutralino mass matrix, and can be 
derived in the following way. By definition of the mixing 
matrix N 1^2]) . we have: 

(TVM)^^. = m,% 
Then for any matrix P, the following identities hold: 



Ofj = {N^2Cw - iViisiy) {cpNj-i + s/jiV^^) + d^i) (17) 



(to; mj) (NPN-'^)^^ = (N {MP + PM) N'^)^^ (22) 
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As a particular case, if we take for P the isospin operator 
that flips the first higgsino sign compared to the second 
one: 



P = T-,= 





-as 



/OO 0^ 

0-10 

Voo 1, 



we recover the gauge independence relation l(2T|) . The par- 
ticular form of the right hand side of this relation then 
follows from the special structure of the symmetric neu- 
tralino mass matrix: 



M 



A C 
B 



where A is a 2 x 2 diagonal matrix reflecting the Majo- 
rana nature of gauginos, i? is a 2 x 2 anti-diagonal matrix 
reflecting the Dirac nature of the SU (2)-charged higgsino 
pair, and C ~ mz is a 2 x 2 matrix that vanishes when 
SU{2) is unbroken and with null determinant to other- 
wise ensure masslessness of the photon. These conditions 
make PM + MP off-diagonal, explaining the non-trivial 
vanishing of the left-hand side of l(22|) with mz- 

These remarks stress the central role played by gauge 
invariance in the structure of the mass matrix. The ap- 
pearance of P = T3 is no surprise, as all xxZ coupHngs 
flnd their root in the gauge invariant hhZ: higgsinos can 
couple to the Z boson thanks to their Uy{^) charge which 
dictates a Dirac behaviour. Only the spontaneous break- 
ing of SU{2) and ?7y(1) carried by C can then split this 
degenerate Dirac system into a pair of Majorana particles. 



Aixx^Zh)^ 



*5 — 



X75X 



(24) 



2 cos-^ Q^fi/ 2mz 

( 2m?y — TO? 
X [1 + ^ 2^ 

The flrst "contact" term in the parenthesis gives a con- 
tribution A ~ -Js = which is divergent and violates 
unitarity in the high energy limit, while the second term is 
better behaved thanks to the appearance of a propagator 
denominator. From a diagrammatic point of view, this is 
expressed by splitting the diagram into a "pinched" part 
and a rest: 

X h X h 

G / 

Rest 




X Z X z 

For the neutralino exchange channel, we have 



A{xx^Zh)^^~2iJ2 



(32 _ jyfi 



xx(^'')(^ + "^xJ7^75X(p) — 

mz 

Expressing ^ = — — to^) — — m^J + (m^ — to^J, 
the flrst term vanishes on-shell, while the second cancels 
the propagator pole to give a contact term and a rest: 



4 High Energy Unitarity 

Despite vague similarities, the gauge independence rela- 
tions l(2T|) do not yet explain the cancellation of ([2|) and 
((H, and in particular their different powers of mixings 
Nij. To get further, we need another relation. Using the 
pinch technique [23], we therefore turn to the high energy 
behaviour of the amplitude. Indeed, from rotation invari- 
ance, the outgoing Z boson must have a pure longitudinal 
polarization, and it is well known that this can lead to 
conflicts with the perturbative unitarity constraint that 
scattering amplitudes be bounded by a constant A(s 
00) < K . 

Having checked gauge independence, we can for sim- 
plicity use the Feynman gauge ^ = 1, to have a purely 
transverse Z propagator free from high energy divergences. 
The dangerous diagrams which must cancel are then the 
s-channel Goldstone and t-channel neutralinos exchanges. 

A Hght-Hke vector being orthogonal to itself, the po- 
larization vector at high energies is approximately 



mz ' 



(23) 



Using this and the kinematic identity 2q.k = q^A-m\ 
the amplitude for Goldstone exchange (fT9l) in this gauge 
becomes: 



A{Zh\ 



-«5 



E 



qZ Qh 



2c^ ^ mz 



X {xibX 



(25) 



Q2 _ j^2^ 

which can be diagrammatically represented by: 
X h X h 



X 




Rest' 



X 



X 



z 



Cancellation of the contact terms in l(24l) and ([Ml 
requires the following identity: 



or, using ([T]), ([4]) and (flTl) . the following relation among 
mixings N: 

(26) 



{N12CW - Niisw) {spNu + cpNi'i) 
4 



1=1 



{{cwNi2 - swNii) [NiiSp - Nacp) 
+ {cwNa ~ swNii) {NiiSp - A^isC/j)} 
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As complicated as it may seen, this equation simply fol- 
lows from the orthogonality condition NijNkj = Sik, with 
i, fc = 3, 4. We have therefore shown that high energy per- 
turbative unitarity of the amplitude is guaranteed by the 
unitarity of the mixing matrix. 

By combining l(26|l and ((20|) , a new non-trivial identity 
among coupHngs is obtained: 



my 



mz 



which is even less trivial when using the definitions of Of 



and O^f 



""x 
mz 



Nt, 



y.{{cwNi2 - swNii) {N^iSp 



(27) 



This identity relates the couplings appearing in Z-exchange 
|[2|) and x-exchange |[3]) , and suggests to rewrite the second 
as: 



and a perturbation 



W 



mz 







~SwC[3 CwCf} 

SwSf3 —CwSf3 



-SwCf3 SwSf} 
CwCp —CwSj} 






triggered by EW-symmetry breaking. Following standard 
Rayleigh-Schroedinger perturbation expansion we start by 
solving the unperturbed eigensystem 

whose eigenvalues are m° = diag (Afi, A/2, —A*, y) and whose 
eigenvectors form the mixing matrix 



/lO 
1 



The first order corrections to eigenvalues m\ — (1^9" | jiys" ) 
are simply the diagonal elements of = N'^W N^'^ : 



A{xx-^Zh)=2iV2p, 



mz cos^Ow 

t{ J 



(28) 



where the first term exactly cancels the s-channel Z ex- 
change, and the second is subdominant in the high energy 
limit s ^ m^,m^. It is however not clear why those re- 
mainders should be subdominant at rest, because even for 
large s = 4m^, there is confusion between dynamically 
suppressed contributions ' 
mixing suppressed ones 
definition 



i^/s <^ 1 and neutralino 
r/m'^ ^ 1. Moreover, the 



„ , (2771^. rriy — 4m^ — 2m?,. 



2m2 



2m^. 



■(29) 



mz 













V2 


V2 








V2 


cws+ 
V2 


V2 


V2 








V2 


V2 









(30) 



where s± = spztcp. From the structure of the perturbation 
W, these diagonal elements clearly vanish. Furthermore, 
first order corrections to eigenvectors 



E 



3'<- 



m^ — m^- 



can be regrouped into N^'^ = (cpj, 'f2Tvh vl) with: 



implies that for mz <C m^, « —R4 are comparable — mz 



with the first term in ([281) . and only their sum becomes 
negligible. To analyze this final cancellation, a perturba- 
tive expansion in mz/m^ is therefore needed. 



5 Perturbation theory 

In the SUSY-decoupling Hmit {mz ^ Mi, M2, fi), the 
neutralino mass in Appendix [C] is naturally split [24] into 
M = Mo + W, with a leading contribution 














V 



%/2(^-Mi) \/2(M2-Ai) 



_Cw_C2_ 






_ _£w_S2_ 









and Ci,2 = (Ats/3 + Mi^2Cp), Si,2 = {Mi^2Sf3 + fJ-cp). 

For a bino-Hke neutrahno {Mi < M2, fi), the Zxx ver- 
tex does not exist without electroweak symmetry breaking 
and knowing the diagonalization matrix N = N° + N-^ up 
to first order in mz allows in fact to compute the first 
non-trivial contribution to Z-exchange ^ which appears 
at second order: 



Mo 



'Ml 
A '/2 

0-^0 



A{xx^ Zh)z^-V2(3, 



Ml 
m% 



'^C2f3t'^ 



{Ml - ^^2) 



(31) 



At the same order, the first non-trivial contribution to %- 
exchange ([3]) requires only one perturbation of the ZxXi 
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vertex, and no perturbation of the hxXi vertex, which does 
exist in the absence of electroweak symmetry breaking. 
Further expanding the propagator in powers of mz would 
give the same structure as ((28|) : 

AP 

A{x^ Zh)^ = V2f3zh^g'^cos2ptan^9w (32) 

rn^_ f l ml + ml \ 
{Ml - /i2) ^ 2 M2 + y ' 

suggesting that the remainder is 0{m%), i.e. two orders 
lower than the leading term. However, to firmly establish 
this conclusion requires a justification for the absence of 
terms 0(m|), which we shall now find by examining the 
general structure of the perturbative expansion. 
When solving the eigensystem 

(Mo + W)^^ Nji = Nlmi (33) 

by power expansions in mz for eigenvalues and eigenvec- 
tors: 

m, = nil +m\ + ... (34) 
N ^ + + ... (35) 

we can to all orders choose the correction to an eigenvector 
in — orthogonal to the corresponding unperturbed 
vector in iV°: the only price is to end up with non-unit 
vectors in N . This choice however simplifies the recurrence 
relation for the solution at order q to: 

m1= {N'^WN'i-^^)^^ (36) 
(N^^m^) = ^ . ^ -YmP n 

V Jji ^0 _ ^0 ^ 1 ^0 _ ^0 

^ J p—1 ^ J 

The expressions for g = 0, and q = I were given above. 
We saw that Mq is 2 x 2 block-diagonal, and so is TV", 
whereas W and thus are block off-diagonal. Following 
the recurrence, this can be generalized, to show that N'^ 
must be block diagonal for q even and block off-diagonal 
for q odd. Because of this structure, mf will vanish for 
q odd, so that mi{mz) is holomorphic in m|. In a sim- 
ilar way, diagonal blocks of TV have a purely even power 
expansion in mz, whereas off-diagonal ones only contain 
odd powers. 

These results can be extended to the amplitudes Az 
and which then contain only even powers of m^: the 
lowest order 0{m^^) vanishes for both as it should to 
allow for a smooth mz limit, the next 0{m'^) is 
equal and opposite for s- and i-channel exchanges, and 
the remainder 0{m\) dictates the amplitude of the Zh 
annihilation channel at rest to be lower than the ti pair 
channel. This can be loosely expressed as: 



A {Zh) Gc—^^A {tt) cx — < A {Zh)y cx 1 

m^ 



The order of magnitude and the power of the suppression 

a{xx^Zh)^ Mf 
then agree with those displayed in figs. I2|3I 



6 Conclusion 

In this work, we have given a quantitative understand- 
ing of why neutralinos at rest cannot annihilate predomi- 
nantly into Zh. The possible relevance of this process for 
indirect DM detection has been shown in the Introduction, 
before pointing out similarities and differences between 
the tt and Zh annihilation channels. We also stressed how 
much a naive estimate of this last channel can fail by ignor- 
ing the subtle but tight Hnks between couplings imposed 
by symmetries, especially broken ones. Because of these 
links, a fine cancellation does occur which requires a closer 
look. Having noticed that this cancellation was getting 
finer with increasingly broken supersymmetry, we showed 
that broken gauge symmetry had to be investigated. This 
was done both at the level of gauge independence in 
gauges, and of unitarity at high energy, known to be deli- 
cate for longitudinal gauge bosons. Both constraints led to 
non-trivial relations among couplings which showed that 
indeed, the SM particles exchanges can be cancelled by su- 
perpartners exchanges, as heavy as these might be. How- 
ever, to quantitatively estimate the importance of what 
remains after this cancellation required to show that a 
perturbative expansion of the amplitudes in mz / m^ con- 
tained only even powers. Whether such cancellation can 
be extended from large s to large m^ = y/s/2 at rest for 
all diagrams suffering from large s unitarity problems, re- 
mains an open question. 

It is a pleasure to acknowledge fruitful discussions with A. 
Djouadi, P. Gay, R. Grimm, J.-L. Kneur, Y. Mambrini, V. 
Morenas, G. Moultaka and J. Papavassiliou. This work was 
partially funded by the GdR 2305 "Supersymetrie" of the French 
CNRS. 



A Indirect Detection Energy Spectra 

The neutrino and photon differential energy spectra of 
figure 1 were extracted from a PYTHIA simulation of 
10^ events for each channel, shown as a function of x ~ 
E^^ -y/m^^. For the hard (anti-)neutrinos from the Z-boson 
in the Zh channel that concerned us most, we have been 
careful to correct the unpolarized PYTHIA results by a 
factor cx x{l — x) translating the purely longitudinal po- 
larization of the Z-boson, which suppresses forward neu- 
trinos with respect to the WW channel. In spite of this 
factor, the Zh channel produces the next-to-hardest neu- 
trino spectrum. 

For photons, the hard (fiat) component of the Zh chan- 
nel around x w 1 comes from h loop-decaying into two 
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photons and therefore only appears at the largest values 
of X. This leaves only a small number of events (~10/bin) 
and large statistical errors which do not appear in the 
fit shown in the bottom figure 1. There are of course no 
events for x > 1, but the precise shape of this vanishing 
(probably similar to that of neutrinos from WW above, 
as shown) is hidden by these errors. 



M = 



Ml 
M2 

mz X C 



mz X C 

,T -/.J 

/i 



(41) 



C is the 2x2 electroweak-breaking contribution to this 
neutralino mass matrix : 



B Dirac and Majorana fermion conventions 

To represent the Clifford algebra of Dirac matrices 
we use the chiral basis: 



cr^ 

CT^ 



75 = 



-1 
1 



where 



cr" = (l,cr) ct'^ = (1,-cr) 

are 4-d extensions of the Pauli matrices cr = (tr^, cr^, cr'^) . 

For Majorana spinors, we use [25] the usual Dirac Feyn- 
man rules after having chosen an arbitrary orientation of 
the spinor lines. The Majorana condition is 



(37) 



and the plane wave expansion of a Majorana field operator 
is thus [26]: 



^^^Y.{-s{p)uAp)e-'^- (38) 



Implementing the Majorana condition ((37| in ((38|) . we find 
the relations 



(39) 
(40) 



which allow to fiip the orientation of external Majorana 
lines. 



C Neutralino Mass Matrix 

The neutralino mass eigenstates are linear combinations 
of gaugino and higgsino fields {ji, W^, Hti Ht^ 

which diagonaHze the mass matrix: 



C 



~SwCf3 swsp 
CwCf3 —Cy/Sp 



with 



sw = smt^^-, cw = cos f)w 
Sf} — sin/3, Cf} — cos/3 

The normalized eigenvectors can be collected into a 
unitary matrix N satisfying: 



N*MN-^ = Mr 



(42) 



where Mu is a diagonal matrix containing neutralino masses. 

In the absence of CP violating phases, the matrix N 
can be chosen as a real matrix, and at least one neutralino 
mass is then negative. 



D Lagrangian terms for relevant couplings 
D.l Z - X - X coupling 



with 



A cos 6 



qA 

4 cos 



w 

9 Ir^Z 



(43) 
(44) 



and 



In the absence of CP violation, there is a basis such that 

cY, = 0. 



D.2 h-Z-Z coupling 



L = -CzhZ^Z^^ 



where 



Cz 



gmz . , , 
■ sm [a — p) 



(45) 



(46) 



COSiJw 

In practice, we will always work in the decoupling limit 
niA ^ fT^z, which implies: 

TT 



a = (3 

Hence, we have sin (a — /3) ~ — 1 in (JH 



(47) 
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D.Z h- X - X coupling 



{K%Pl + K.jPr) } xj 



with 



L, 



Na {Nj2 - taii9wNji) + (^ ^ j) 



The symmetric form of K and L comes from the fact 
we are working with Majorana particles: 



Xi (1 ± 75) Xj = Xj (1 ± 75) Xi 
For real N, this Lagrangian simplifies to: 



^ = J E 



Xj 



where 



and 



2 cos 9w 



O'lj = {cwN,2 - swNa) (saTV^s + c^^Nji) + (i <-> j) 
D.4 Z-fermion-fermion coupling 



where 



r<ZA 



2 cos dw 
9 



2 cos 6* 



w 



r3/-2sin2 9wQf) 



Q f and T3 / are the charge and third weak isospin compo- 
nent, with the usual normalization: T^top = 1, and Qtop — 
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